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Abstract 

We analyze the CP problem in the context of a supersymmetric extension of the stan- 
dard model with universal strength of Yukawa couplings. A salient feature of these models 
is that the CP phases are constrained to be very small by the hierarchy of the quark masses 
and the pattern of CKM mixing angles. This leads to a small amount of CP violation 
from the usual KM mechanism and a significant contribution from supersymmetry is re- 
quired. Due to the large generation mixing in some of the supersymmetric interactions, 
the electric dipole moments impose severe constraints on the parameter space, forcing the 
trilinear couplings to be factorizable in matrix form. We find that the LL mass insertions 
give the dominant gluino contribution to saturate Ek- The chargino contributions to e'/s 
are significant and can accommodate the experimental results. In this framework, the 
standard model gives a negligible contribution to the CP asymmetry in B-meson decay, 
aj/ipKs- However, due to supersymmetric contributions to B d — B d mixing, the recent 
large value of aj/^K s can be accommodated. 
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1 Introduction 



The understanding of the origin of fermion families and the observed pattern of fermion 
masses and mixings, together with the origin of CP violation, are among the major out- 
standing problems in particle physics. In the standard model (SM), since the flavour 
structure of the Yukawa couplings is not defined by the gauge symmetry, fermion masses 
and mixings are arbitrary parameters, to be fixed by experiment. A salient feature of the 
pattern of quark masses and mixings is the fact that the spectrum is strongly hierarchical 
and the mixings are small. 

Concerning CP violation, it arises in the SM from complex Yukawa couplings which 
lead to a physical £>ckm phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The 
strength of CP violation in the SM is given by Im |Qc.p|> with Qcp denoting any rephasing 
invariant quartet of the CKM matrix. If one considers the quartet Qcp = (KsKbKTfcKs); 
the experimental value of the moduli entering in this quartet constrains |Qcp| to be 
of order 10~ 5 , while the observed strength of CP violation in the kaon sector requires 
(f) = aig{y us V c bV* b V* s ) to be of order one. Although in the SM no explanation is provided 
for the experimentally observed values of |V^|, once these values are incorporated, the 
SM accommodates in a natural way the observed strength of CP violation in the kaon 
sector, measured by the parameter ek- On the other hand, it has been established P 
that the strength of the CP violation in the SM is not sufficient to generate the observed 
size of the baryon asymmetry of the Universe. This provides the strongest motivation to 
consider new sources of CP violation, beyond the one present in the SM. 

In supersymmetric (SUSY) extensions of the SM there are additional sources of CP 
violation, due to the presence of new CP violating phases. However, these new phases are 
severely constrained by experimental bounds on electric dipole moments (EDMs) IS] ■ 
Setting these new phases to zero is not natural in the sense of 't Hooft jlj since the 
Lagrangian does not acquire any new symmetry in the limit where these new phases 
vanish. This is the so-called SUSY CP-problem. Since the question of CP-violation is 
closely related to the general flavour problem, one may wonder whether it is possible, 
within a supersymmetric extension of the SM, to establish a connection between the need 
for small CP violating phases and the observed pattern of quark masses and mixings. 

Such a connection is possible if one imposes universality of strength for Yukawa cou- 
plings (USY) on a SUSY extension of the SM. The idea of USY consists of assuming that 
all Yukawa couplings have the same modulus, the flavour dependence being all contained 
in the phases [5J . This leads to what is known as pure phase mass matrices [B] . This class 
of couplings can be motivated by local horizontal symmetries (see Appendix), and it has 
also been recently shown that such textures may naturally arise within the framework of 
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a theory with two "large" extra dimensions jjj- 

The fact that the quark masses are strongly hierarchical, together with the smallness 
of \Vcb\, \V u b\ implies that these phases have to be small, at most of order 1712/7113 [13 IB]- 
One is then led to a scenario where the strength of CP violation through the Kobayashi- 
Maskawa (KM) mechanism is naturally small, in fact too small to account for the observed 
CP violation (CPV) in the kaon sector. This provides motivation for embedding the USY 
ansatz in a larger framework, where new sources of CP violation naturally arise. 

In this paper, we impose the USY ansatz in the framework of SUSY. This extension is 
specially interesting since new SUSY contributions to CP violation in the kaon sector may 
play an important role in providing sufficient CP violation to saturate the experimental 
value of Sk- 

We study the question of electric dipole moments (EDMs) in supersymmetric USY 
models. Although the CP violating phases of the Yukawa couplings are small, the USY 
ansatz implies a large generation mixing in some of the SUSY interactions, which in turn 
leads to potentially large contributions to the EDMs. The recent bound from the mercury 
atom EDM jHj sets important restrictions on the class of models we are considering. We 
will show that due to these constraints one can only saturate the observed value of Ek 
with non-universality of the soft trilinear couplings. 

Finally, we analyze CP violation in the B-sector. Both the BaBar ^Hj and Belle jTTj 
collaborations have provided clear evidence for CP violation in the B-system, although at 
present the experimental errors are relatively large. The large CP asymmetries observed 
by BaBar and Belle are in rough agreement with the SM predictions. However, there is 
still quite a large room for new physics and models with small CP violating phases like 
USY, as well as some of the models with CP as an approximate symmetry jT2j might 
remain phenomenologically viable, provided there are new physics contributions to Ek, 
e'/e, and B d — B d mixing [TB"] . 

This paper is organized as follows. In Section 2 we present the main features of the 
USY ansatz for the Yukawa couplings. We show that USY is preserved under the renor- 
malization group (RG) evolution, and give a specific example of a choice of USY phases 
that leads to a correct quark mass spectrum and |V^ KM |. In Section 3 we implement 
the USY hypothesis within a supersymmetric extension of the standard model and in- 
vestigate how the EDM constrains the SUSY parameter space. The remaining of the 
section is devoted to the study of the kaon system CPV observables. In Section 4 we 
analyze the implications of small CPV phases for CP asymmetries in Bd, meson decays. 
Our conclusions are presented in Section 5. 
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2 Universal strength of Yukawa couplings 



In the absence of a fundamental theory of flavour, one is tempted to consider some special 
pattern for the Yukawa couplings, in the hope that they could naturally arise from a 
flavour symmetry, imposed at a higher energy scale. As mentioned in the introduction, 
one of the simplest suggestions for the structure of the Yukawa couplings consists of 
assuming that they have a universal strength g y , so that the mass matrices can be written 
as rriij = g y v / V2 exp(i0jj) , where i,j denote family indices. In the limit of small phases, 
the USY matrices can be viewed as a perturbation of the democratic type matrices [Hj. 

Within the USY framework, the quark Yukawa couplings can be parametrised as 
follows: 

Uu = exp and D iA = -4 exp , (1) 



L -- lJS — tJ g x- j 

where X U) d are overall real constants, and $ M,d are pure phase matrices. By performing 
appropriate weak-basis (WB) transformations the matrices U and D can, without loss of 
generality, be put in the form 
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V 



1 \ 



D 



J 
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1 
1 



1 \ 



K 



(2) 
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where K = diag(l, e tK1 , e lK2 ). 

The phases «j do not affect the spectrum of quark masses, but influence Vckm- The 
phases (p,q,r,t) affect both the spectrum and Vckm- As previously remarked, the hier- 
archical character of the quark masses constrains each of these latter phases to be small. 
In order to see how these constraints arise, it is useful to define the Hermitian matrices 
H u ,d = 3/A^ d QQ\ where Q = U,D. It can then be readily verified that the second 
invariant of H is given by j3] 



X(H) 



4^ . 

— / si 
9 ^ 



sm 2 (T r 
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m\m\ + m\m\ + m|mg 
(mf + m\ + m\) 2 



(3) 



where, in each sector, a m are linear combinations of the phases (p,q,r,t). Since x{H) is 
a sum of positive definite quantities, the modulus of each of the phases has to be small, 
at most of order 777.2/ '^3- As mentioned above, the phases Ki do not affect the quark 
spectrum. However, it has been shown jH] that the fact that (|V c fe| 2 + | V^b] 2 ) is small, 
constrains each one of the phases Ki, k 2 to be at most of order |V^,|. As we will see, these 
k phases will prove to be essential to control the size of the Vckm elements, in particular 
that of V t d- 
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In our search for a choice of phases leading to the correct masses and mixings, we 
expand both x{H) and the determinant of H in the limit of small phases, and per each 
of the quark sectors, we write two of the phases as functions of the masses and of the 
remaining two phases. At this point, we have still some freedom, which we use to fit the 
experimental bounds on the Vckm matrix. Therefore, for fixed values of the quark masses* 
and using as input the phases p u ,d, fu,d and Ki, we find regions of the phase-parameter 
space where one correctly reproduces the masses and mixings of the quarks. 

As referred in the introduction, USY can be generated at some high energy scale 
through the breaking of a flavour symmetry. In Ref. 0, it was shown that USY can 
account for the quark mass spectrum at low energy scales and thus it is important to 
investigate whether this structure for the Yukawa couplings is preserved by the renormal- 
ization group equations (RGE) [16 . 

At the weak scale, the phases required to explain the fermion masses and mixings are 
small, hence we can expand the Yukawa couplings of Eq. (JTJ) as 



X 



f [E tJ + <$« ] and Dij = ^ [E tJ + . 



(4) 



In the above equation, we neglected terms 0(Q 2 ), and E is the democratic matrix, defined 
as Eij = 1, Vz,j . The RGE for the real part of Eq. (jlj) can be easily identified with the 
evolution of the top and bottom Yukawa couplings: 

d X u 1 



dt 

d X d 



16vr 2 
1 



(-G u + QX 2 u + X 2 d ) X u , 
(-G d + QX 2 d + X 2 u + X 2 ) X 



dt 16tt 2 v ~" ' ~" d ' " u ' " t/ ""' (5) 

where Gkit) = Y^=i c l 9i(t)> with gi(t) the running gauge couplings for SU(3) x SU(2) x 
U(l) gauge group and Ck the MSSM RGE coefficients. The reason why, in this limit, 
the running of the overall Yukawa coefficient decouples from the running of the phases 
becomes clear if, by means of a weak basis transformation, one moves from the USY basis 
to the "heavy" basis (Y Y' = F^YF), where 

A 



Y' 



diag(0,0,3)+^ 



u,d 



with 



(6) 



^ k (t) = F^ k F and F 



1/V2 -l/y/6 1/V3 
2/y/G l/y/3 

-l/y/2 -1/V6 1/V3 



(7) 



*At ji = m t0 p, we use the following values for the running quark masses: m u = 0.0018 GeV, nid — 
0.0048 GeV, m s = 0.075 GeV, m c = 0.6 GeV, m b = 2.9 GeV, m t = 170 GeV, consistent with the values 
given in Ref. |15| . 



The lighter generation Yukawa couplings are thus small perturbations of the leading third 
generation, originated by the small USY phases. 
The running of the imaginary parts is 



d 

dt 
d $ d 
dt 



16vr 2 
1 

16tt 2 



A d, 



Xi(-3I + E) + -^(-37 + E) 
X 2 d {-31 + E) + ^(-37 + E) 



167T 2 
1 

16tt 2 



E 



\2 \2 
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lM A d E + —A U E 
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where J denotes the 3x3 identity matrix and A k = & k E — E(§ k ) T . The scale dependence 
of the phases can be explicitly analyzed if we write Eqs. (jHJ) in the "heavy" basis: 
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dt 
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(9) 



(10) 



where we defined k 2 (t) = A 2 (t)/167r 2 . 

We have performed the analysis of the scale dependence of the phases in the framework 
of the SM and in the MSSM. The differential equations were solved numerically and 
analytically, the latter in the limit of small phases. In both cases we have found that if 
the USY phases are small at some initial scale (which in a phenomenologically consistent 
model is justified by the smallness of the phases required to fit the quark masses and the 
Vckm at the weak scale), then the approximation of Eq. (J3J) remains valid to energies up 
to the GUT scale. 

The explicit dependence of the phases on the scale can be easily displayed for the 
low tan/? regime, where one can neglect terms in A 2 , and A 2 relatively to the leading A 2 
contributions. In this limit, the evolution of the up phases is given by 



= a(t)$ u (0) + (1 - a(t)) X u , 
where a(t) = exp ^—3 J * k 2 (t')dt'^j contains the scale dependence and 



(11) 
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^ 3 (0) 

^3(0) 
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(12) 
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is a function of the initial phases; the structure of \ u allows these to be combined so that 
the up quark Yukawa coupling for low tan j3 values takes the form: 



U(t) 



K u\ e ia(t)* tt (0) K u 



R 5 



(13) 



where K"£ and K\ are diagonal phase matrices, whose entries depend on the initial phases 
contained in the matrix \ u an d on the integral a(t). 

The analysis for the down sector is entirely analogous. In the low tan/? regime, the 
solution for the running of the down phases takes the form: 



= b{t)¥{0) + (1 - b(t)) X d + d(t)r} u 



(14) 



where b(t) = exp y — f kl(t')dt'j and d(t) 
are also functions of the initial phases. 



t kl(t') 
6(f) 



dt' . The matrices x d an d ^ 



X 







vixtp) ^ 2 (o) ^ 3 (o) 



Ft 



^ 3 (0) 
^ 3 (0) 




Ft . (15) 



Just like in the previous case, these phase matrices can be reorganised so that the down 
Yukawa coupling evolution reads: 

Xd(t) 



D(t) 



Kf Q e ib ®* d W K d R 



(16) 



where Kf, and Q are diagonal phase matrices; the first two depend on the function 
b(t) and on the phases appearing on x d , while the third depends on d(t) and r] u . 

By means of yet another weak basis transformation, one can still rewrite the Yukawa 
couplings in Eqs. (|13I16|) as 



U(t) 



D ia(t)$"(0) . 



D(t) 



X d (t) 



pt ^m d (o) p (p = Q\ K d L K u J }) 



(17) 

This means that it is always possible to move to a basis where the Yukawa couplings 
are easily related to the original ones. Obviously, the P matrices have no effect on the 
quark spectrum, only on Vckm- Therefore RG evolution preserves the USY texture that 
is responsible for the non-degenerate quark spectrum. For the case of large tan/3, the 
argument can be equally verified. 

At this point, one should make a final remark regarding the possibility of having Her- 
mitian and USY Yukawa couplings, that is having the phase matrices verify the condition 
c&lj = — (q = u,d). As it was first pointed out in [T7j . Hermitian USY leads to a 
sum rule for the quark masses, namely mi ni2 + vn^ + mi = 0, which is clearly 



violated by experiment. The question that naturally follows is whether Hermitian USY 
textures imposed at some higher energy (GUT scale, for example) can produce a phe- 
nomenologically viable texture through RG effects. As it can be seen from Eqs. fll3pi6) 
117)1 . if one starts with Hermitian structures, that is $^ (0) = — $^(0), we find that at any 
other scale one can always rotate to a basis where the phases that govern the spectrum 
scale with the flavour-independent coefficients a(t) and b(t), so that the impossibility of 
having Hermitian USY textures holds at any scale. 

To illustrate some of the features of this ansatz, let us consider, as an example, the 
following choice of USY phases: 



U 



D 



K 



~3~ 



( 



,t 0.000517 



-i 0.01438 



A 0.00229 



-i 0.09278 



-i 0.0000074 

1 
1 

A 0.02833 



1 \ 

A 0.00158 



V 



diag(l , 1 , e 



1 

i 0.0138\ 



-i 0.1453 



K 



(18) 



The overall factor has been defined as A„/3 = m t /v sin/3 and Ad/3 = ra^jv cos /?. 
The Yukawa couplings can be diagonalized by means of the transformations 



y; 



d{u) 



diag ' 



(19) 



where Sr and Sl are unitary matrices. In the present USY model we find 



SI 




0.0236- 1.66 x 10 
0.816-3.72 x 10" 



-4 



-0.719 + 1.83 x 10- 4 i \ 
-0.388 + 1.65 x 10~ 4 i 



(20) 



-0.264 
\ 0.578 i 
and the associated Vckm 

/ 0.975- 6.88 x 10~ 4 % 

Vckm = 



2.56 x 10~ 3 + 0.577 i 9.79 x 10" 5 + 0.577 i J 



( 0.772 -0.157 + 0.0122 i -0.615 - 4.99 x 10~ 4 i \ 
0.798 - 0.0759 i -0.537 + 3.05 x 10" 3 i 



(21) 



0.511 + 0.575 i -2.07 x 10~ 3 + 0.577 i J 



\ 0.0063 + 1.10 x 10 



-4 



0.221 + 3.96 x lO" 3 i 

0.973 + 0.061 i 
-0.041 - 1.75 x 10" 3 j 



0.0027-4 x 10" 
0.041 + 1.92 x 10 



-3 



0.999 + 0.027 i 



(22) 

Within the context of the SM, the matrices S^' R do not have, by themselves, any physical 

~~ Vckm is physically meaningful. However the 
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meaning, only the combination S^.S^ 



matrices S^' d R do play a significant role in some extensions of the SM, as, for example 
the MSSM with non-universal soft-breaking terms. We shall discuss this question in the 
following sections. 

The values of the moduli of the above presented Vckm are i n agreement with the 
experimental bounds. The central value of \V t d\ {\Vtd\ ~ 0.0063) is somewhat smaller than 
the central value of \Vtd\ extracted from experiment ^Hj, within the context of the SM, 
namely \Vtd\ = 0.0073. However, it should be noted that in the framework of the SM, 
the value of \Vtd\ is derived from the experimental value of Bd — Bd mixing, and not from 
tree-level decays, as it is the case for the moduli of the first two rows of Vckm- Since 
in the SM, Bd — Bd mixing only receives contributions at one loop level, there may be 
non-negligible contributions to Bd — Bd mixing from physics beyond the SM. 

Concerning CP violation, it can be readily verified that the above Vckm leads to a 
strength of CP breaking (measured by J = Im {V u dV cs V* s V* d ) [T^j) that is too small to 
account for the observed value of Sk- In fact, for the above example, one obtains J ~ 10~ 7 , 
which is to be compared to the required value of J ~ 10~ 5 . This is a generic feature of 
the USY ansatz. Although one can obtain somewhat larger values of J (e.g. J ~ 10 -6 ) 
for a different choice of USY phases j^J, the tendency in this class of ansatze is to have 
values of J which are too small to saturate the experimental value of Ek- This motivates 
the embedding of USY in a larger framework where new sources of CP violation naturally 
arise. In the following sections we address this question in the context of supersymmetric 
extensions of the standard model. 

3 Supersymmetric USY models and CP violation 

As advocated in the introduction, supersymmetric extensions of the SM may provide a 
considerable enhancement to CP violation observables, since in addition to containing 
new "genuine" CPV phases, they introduce new flavour structures. Accordingly, SUSY 
emerges as the natural candidate to solve the problem of CP violation inherent to the USY 
models discussed in the previous section. We recall that within the USY framework the 
amount of CP violation arising from the CKM mechanism is typically quite small, hence 
one will need significant supersymmetric contributions that both succeed in saturating 
the values of the CP violating observables and in having EDM contributions which do 
not violate the current experimental bounds. 

We will consider the minimal supersymmetric standard model (MSSM), where a min- 
imal number of superfields is introduced and R parity is conserved^, with the following 

^For models with broken R parity where universal strength of Yukawa couplings is also assumed, see 
for example Ref. |2(Jj . 
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soft SUSY breaking terms 



V SB = mlMa + e ab (A^H b 2 qlu% + A%Y^ Li d\ + A[ j Y i l j H"l b L .e* R , 

- BfiH^Hz + H.c.) - ^{m 3 gg + m 2 W a W a + m x BB) , (23) 

where i, j are family indices, a, b are SU(2) indices, and e a b is the 2x2 fully antisymmetric 
tensor, with ei 2 = 1. Moreover, <p a denotes all the scalar fields of the theory. Although 
in general the parameters fi, B, A a and m ; can be complex, two of their phases can be 
rotated away |21j . In a minimal SUGRA scenario (mSUGRA), the soft SUSY breaking 
parameters are universal at some very high energy scale (which we take to be the GUT 
scale), and we can write 

ml a = m 2 0} m l = m 1/2 , A% = A e^ . (24) 

In this case, there are only two physical phases 

(j) A = arg(A*mi /2 ), M = arg(^ m 1/2 ). (25) 

In order to have EDM values below the experimental bounds, and without forcing the 
SUSY masses to be unnaturally heavy, the phases 4>a and </> M must be at most of order 
10~ 2 [2 j. Thus, in this class of SUSY models with minimal flavour violation, complex 
Yukawa couplings leading to a physical Sckm phase in the CKM matrix are the main 
source of CP violation. It should be also stressed that even if one ignores the bounds 
from the EDMs, and allows <pA and 0^ to be of order one, the latter models do not 
generate any sizeable new contribution to Ek and e'/e [22] apart from those present in 
the SM. In the mSUGRA scenario of Eq. (j24|) . the USY structure of the Yukawa couplings 
is inherited by the trilinear terms so this simple flavour structure is excluded, since, and 
as aforementioned, no new contributions to CPV observables are generated. However, 
in general supergravity scenarios [23]; the Yukawa couplings can be independent of the 
hidden sector fields that break supersymmetry and hence one can have A-terms with a 
flavour structure that is not related to that of the Yukawa couplings. This has recently 
motivated a growing interest on supersymmetric models with non-universal soft breaking 
terms, and a considerable amount of work has been devoted to the analysis of the effects 
of the new flavour textures on the CPV observables [2H I2SI I2S] • 

Regarding the non-universality of the trilinear soft terms, it has been emphasized that 
a new flavour structure in the A terms can saturate the experimental bounds on ek and 
e'/e even in the presence of a vanishing Sckm 1231123 However, the phases associated 
with the v4-term diagonal elements might induce large EDMs. The latter problem can 
be overcome if one takes the Yukawa couplings and the A-terms to be Hermitian, in 
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order to ensure that the diagonal elements of AijYij are real in any basis j2H 1211 • This 
problem is also less severe if the trilinear terms can be factorized as A^Y^ = A.Y or Y.A, 
an interesting possibility that naturally arises withing the context of general SUGRA 
models [22] • As we will discuss below, this factorization implies that the mass insertion 
(Sfx)LR (which is the most relevant to the EDMs contributions) is suppressed by the ratio 
rrid/mq, where is the down quark mass and rriq is an average squark mass. 

The non-universality of the squark mass matrices, which can be arbitrary in generic 
SUSY models, is not constrained by the EDMs. Nevertheless, Amx and Ek impose 
severe constraints on the squark mixing, which can still be avoided in models with flavour 
symmetries or squark-quark alignment [2*%j . 

After this brief discussion, we begin the analysis of the implications of having USY 
Yukawa couplings within a supersymmetric scenario of CP violation. As in most SUSY 
models, the trilinear terms will play a key role, since they are directly related with the 
increasingly more stringent bounds coming from the EDMs. In what follows, we will 
study the enhancement of the CP observables through the non-universalities of the soft 
breaking terms in supersymmetric USY models. 

3.1 Constraints from the electric dipole moments 

In this section we investigate how the EDM bounds constrain the parameter space for 
SUSY models with universal strength of Yukawa couplings. The current experimental 
bound on the EDM of the neutron is given by [2H] 



This bound can be translated into constraints for the imaginary parts of the flavour 
conserving LR mass insertions. 



is assumed to be real and its magnitude is computed from electroweak symmetry break- 
ing. Recall that and are the rotation matrices that diagonalize the down (up) 
Yukawa couplings. For the specific case of USY couplings discussed in Section El the 
structure of S^ u ^ was presented in Eqs. (j2QI21j) . manifest the same large intergen- 

erational mixing and large phases. The effects of these phases are absent in the SM and 
in supersymmetric models with universal soft-breaking terms (i.e., flavour structures that 
are aligned with the Yukawa couplings). Nevertheless, the effect of any minimal devia- 
tion from this scenario, particularly regarding the non-degeneracy of the trilinear terms, 



d n < 6.3 x KT 26 e cm (90% C.L.). 



(26) 




(27) 



where m 2 q is the average squark mass and u» 



(Hf) /a/2. In our analysis, the \i term 
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is strongly enhanced due to the large mixing, providing significant contributions to CP 
violation observables. On the other hand, the bounds arising from the experimental mea- 
surements of the EDMs severely constrain any non-degeneracy of the trilinear terms, even 
in the limit of vanishing SUSY phases. Imposing that d n does not exceed the experimental 
limit requires Im (5f[ )lr to be less than 10~ 6 . However, a stronger constraint on these 
mass insertions comes from the recently measured electric dipole moment of the mercury 
atom |Hj 

d Hg < 2.1 x 1(T 28 e cm . (28) 

This constraint corresponds to having Im(5fj )lr ^ 10~ 7 — 10~ 8 . Since the mercury 
EDM also receives considerable contributions from strange quarks, one has in addition 
lm(5 22 )LR < 10~ 5 — 10~ 6 j2]. These bounds impose a stringent constraint on most super- 
symmetric models, implying that the CP violating phases are very small and (or) that 
SUSY soft breaking terms must have a special flavour structure. 

Small CP phases can be motivated by an approximate CP symmetry and as previously 
pointed out, USY provides a natural scenario where all Yukawa coupling phases are bound 
to be < 1CT 1 . In this work we take as a guideline the assumption that all supersymmetric 
phases should be no greater than the largest of the phases present in the Yukawa textures 
of Eqs. (jTgj). 

Although in many models the trilinear terms do not have the same structure of the 
Yukawa couplings, when introducing the USY ansatz in a supersymmetric framework 
the most appealing scenario would be to have a USY texture in the A-terms as well. 
Nevertheless, we have verified that such textures for the A-terms do not suceed in pro- 
viding sizable SUSY contributions to the CPV observables without violating the EDM 
bounds. Therefore we consider a more general scenario for the trilinear terms, and allow 
for non- universalities in the magnitude of the several elements. We will show that, even 
in the limit of small (or vanishing) supersymmetric phases, the large mixing inherent to 
USY couplings (displayed in Sl and Sr) , together with the experimental limits on EDMs, 
severely constrain any non-universality of the trilinear terms. In Fig. we present the 
constraints from the EDMs on the off-diagonal entries of the A-terms, in particular A\ 2 
and Ai 3 . In our analysis we assumed tan/3 = 5, m = mi/ 2 = 250 GeV and = ttiq for 
all elements except A 12j i3, which we set in the range [— 3m ,3m ]. 

From these figures, it is clear that any significant non-universality among the A terms 
(e.g. at 10%) leads to unacceptably large EDMs. Similar constraints hold for the other 
off-diagonal elements A 2 \, v4 31 , A 23 and A 32 . It is also worth noticing that such constraints 
are far more severe than those obtained in the case of hierarchical Yukawa couplings, as 
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Figure 1: Neutron and mercury EDMs as function of the off-diagonal entries A\2 and A13 
for tan/3 = 5 and m = mi/2 = 250 GeV. Aij = 1 Vzj 7^ 12(13) on the left (right) figure. 



a consequence of the already mentioned large intergenerational mixing. 

As recently shown, this problem can be overcome if one takes the Yukawa couplings 
and the A-terms to be Hermitian, in order to ensure that the diagonal elements of AijYij 
are real in any basis [2E]. Another interesting possibility is to have trilinear terms that 
can be factorized as A = A^Y^ = A.Y or Y.A [22] • As above referred, this factorization 
implies that the mass insertion (Sf^LR is suppressed by the ratio mj/m^. 

As an example for the trilinear couplings that avoids the EDM constraints and implies 
signficant SUSY CP violation, we consider the following structure: 



^ a a a ^ 
A = m b b b 

\ c c c J 

In this case the trilinear couplings A can be written as 



(29) 



A = diag(a, b, c) ■ Y . (30) 

With this structure for the trilinear couplings, the EDMs do not impose any constraint 
on the parameters a, b, c if these are real, and the associated EDM values are similar to 
those obtained with universal and real A terms. For complex a, b, and c, the constraints 
on the phases ip a , <fb, and ip c are shown in Fig. 2. 

As it can be observed for the case of universal \Aij\, compatibility with experiment 
requires 

^ a <0.05rad Vb < tt/2 ; y? c < 0.07 rad . (31) 

13 




Figure 2: The mercury EDM as function of the trilinear SUSY phases tp for tan/? = 5, 
tjiq = mi/2 = 250 GeV. \Aij\ = 1 



The reason why the phase ipj, is much less constrained than ip a and ip c is that the mixing 
between the first and the second generations, as observed from Eqs. (I2()|21jl . is smaller 
than the mixing between the first and the third generations. For \A^j\ ~ 3mo, the EDMs 
constraint on these phases becomes more stringent: 

y? a <0.02rad , (p b < 0.35 rad , y? c < 0.02 rad . (32) 

Thus, we can conclude that in USY models the maximal allowed non-universality for 
the A-terms is the structure presented in Eq. fj2Uj) . with the associated SUSY phases (ipij) 
within the limits given in Eq. (|32jl . In view of this, we analyze the contributions to the 
kaon system CP V observables, £% and e' / e. 



3.2 SUSY contributions to ek and e'/e 

We proceed to investigate how CP violation can further constrain the parameter space 
for the case of supersymmetric USY models. In the kaon system, the most relevant 
observables are Ek arid e'/e. We first turn our attention to the indirect CP violation 
parameter, Ek, which is given by 

E ^7fe^ ImMl2 ' (33) 

where Myi is the off-diagonal entry in the kaon mass matrix and ArriK is the mass dif- 
ference of the short and long-lived kaon states. The kaon mass matrix element can be 
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defined as follows 

2m K M* 12 = (K \H* S=2 \K°) , (34) 

where Ti^p 2 is the effective Hamiltonian (EH) for AS = 2 transitions, which can be 
expressed in the operator product expansion as 



njAS=2 
H eff 



x;ci(/*)& 



(35) 



In the above formula, Cj(^t) are the Wilson coefficients and Qi are the EH local operators. 
The relevant operators for the gluino contribution are [SUj 



Qi 
Q 4 



Q 3 



u R b L R Li 



a R b L u L b Ri 



Q§ — ^R S L^X S R> 



(36) 



as well as the operators Qi,2,3, that are obtained from Qi^,3 by the exchange L <-> R. In 
the latter equations, a and f3 are SU(3) C colour indices, and the colour matrices obey the 
normalization Tr (t a t b ) = 5 ab /2. Due to the gaugino dominance in the chargino-squark 
loop, the most significant x ± contribution is associated with the operator Qi |31j . 

In the presence of supersymmetric (g and x^) contributions, the following result for 
the amplitude M 12 is obtained: 



M 



12 



Mf 9 M + ML + M* 



[ 12 



12 



'12 



The SM contribution can be written as [321 EB] 



m sm = ^kf KmK B K r 



(37) 



(38) 



where fx is kaon decay constant, Bk is the bag parameter, and the function T can be 
found in [32]. For the specific USY ansatz we are considering, we find that the SM 
contribution to Ek is O(10~ 5 ). The supersymmetric term Mf 2 is given by 



Mf 2 



-as 1 
216m,|2 3 



12) xx + ^! (a*) (24x/ 6 (x) + 66/ 6 (x 



+ (^l 2 )xx(^ 2 )iiH 



m s (//) + m d (/z) 



5 4 (/i) 378o;/ 6 (x)-54/ 6 (a;; 



+ B 5 (ji) [6xf 6 {x) + 30f 6 (x] 



m K 



255 9 \ 

--BM - -BM) xh { X ) + (5f 2 ) LR (S d u)RL Km ^ )+mM 



-99B 4 (/i) - 455 5 (/x) 



m s (/i) +m d (/j,) 
m K 




(39) 
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where x = (rrig/mq) 2 and the functions fe(x), fe(x) are given in Ref.[30j. We have used 
the matrix elements of the normalized operators Qi(/J>) with the following ^-parameters 
at 11 = 2GeV j3D] 

B 1 (fi) = 0.60 , B 2 (fx) = 0.66 , B 3 (/i) = 1.05 , 

Bi(ji) = 1.03 , B 5 (p) = 0.73 . 

Finally Mjf is HD] 

M t = t^A \m K f K BM (j2 K *aM L ) ab K bl ) Wim^H^x,), (40) 

7bSir z m S 2 3 \ *—f j 

H \ a,b / 1,0 

where Xi = (m^f/rriq) 2 . Here K is the CKM matrix, a,b are flavour indices, i,j label 
the chargino mass eigenstates, and V is one of the matrices used for diagonalizing the 
chargino mass matrix. The loop function H(xi,Xj) is given in Ref. 31 . 
To saturate Ek from gluino contributions, one should have 

^M*? 2 )£ L |~10- 3 or ^ImO^I-lO- 4 , (41) 

whilst the chargino contributions require 

Im(<%)L ~ 10- 1 . (42) 

In the present scenario, using the USY texture of Eq. (JTSj) and A-terms as in Eq. fl2T?j) . 
we find the correlation between these two mass insertions is as given in Fig. El 



From Figure El one can see that the LL mass insertions give the dominant gluino 
contribution to saturate contrary to the LR mass insertions, which are typically 
much smaller than the required value given in Eq. ()41jl . In addition, we also find (for 
m = m\/2 = 250 GeV) that Im(5" 2 )^ ~ 10~ 3 , two orders of magnitude below the 
required value, so that one has a negligible chargino contribution to Ek- Therefore, 
with USY Yukawa couplings, the MSSM with non-universal A-terms (as those given in 
Eq. (|29|l) can account for the experimentally observed indirect CP violation in the kaon 
system. 

It is important to note that the behaviour of the LL mass insertion (particularly the 
significant dependence on the non-universality of the A-terms) is in clear contrast with the 
case of standard hierarchical Yukawa couplings. It is well known that in these scenarios, 
and with universal scalar soft masses, lm(5f 2 )LL is typically of order 10~ 5 and the LL 
mass insertions are independent of the A-terms. In our model, the large mixing displayed 
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Figure 3: Correlation between (Sf 2 )LL an d ($12) lr mass inserions for for tan/3 = 5, 
tuq = mi/2 = 250 GeV. a, 6, c vary from [—3, 3] and (p a ,b,c are chosen so that the EDMs 
are consistent with the experimental limits. 



by the rotation matrices enhances any non-diagonal contribution to M~ induced by 

W L 

the trilinear terms from RG evolution, hence LL mass insertions can account for the 
experimental value of Sk- 

Next we consider SUSY contributions to e'/e. The grand average of the experimental 
values is EE] , 

Re (e'/e) = (19.2 ± 2.4) x 10" 4 . (43) 
e' can be given as ratio of isospin amplitudes Aj, which are defined as 

A ie ^ = (n7r(I)\Ht f S f =1 \K ) , (44) 

where / = 0, 2 is the isospin of the final two-pion state. Therefore we can write 

£ 1 „ _L e ^s ) Im (jg-jo) (45) 
e V2 \eA \ 2 

Assuming that the strong phase difference (S 2 — Sq) is close to 7r/4, the expression for 
e' can be rewritten as f 

L = ji-^—£(n-i) (46) 

e V2\e\ 

with 

uj = ReA 2 /ReA , £ = and fl= . (47) 

Re v4 culm A 
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The real parts of both amplitudes are experimentally known to be Re A% = 1.5 x 10 GeV 
and ReA = 33.3 x 1(T 8 GeV. 

The effective Hamiltonian for AS = 1 transitions receives contributions from operators 
with dimension four, five and six, and can be symbolically written as 

K A / =1 = + «2 + «2 ■ ( 4 8) 
The AS 1 = 1 effective Hamiltonian containing dimension-6 operators is given by |37| 

10 

W$ =^C i (/i)0, + H.c. , (49) 
i=i 

where Cj are the Wilson coefficients. refers to the current-current operators, to 
the QCD penguin operators and O^^iq to the electroweak penguin operators. In addition, 
one should take into account Ot operators, which are simply related to Oi by the exchange 
L <-> R. The corresponding matrix elements (inr\Oi\K ) can be obtained from the matrix 
elements of (7m\Oi\K°) , multiplying them by (—1) while the Ci are obtained from Cj by 
exchange L <->• R. In principle, there are also two dimension-5 "magnetic" operators 7 
and O g [SB], which are induced by gluino exchanges: 

= C 7 7 + C g O g + H.c. , (50) 

where 



a 



a 



QV7T 

7 2 

9 9 



(5f 2 ) iL ^M 3 (x OT ) + {S&LR-f- -Mi(x gq ) 



K 2 ) LL f-~M 3 (z 93 ) - 3M 4 ( % )j + (S&lr 7 ^- (-\M.ixg,) - 3M 2 ( % )j 



(51) 



Here x gq = rn^/m 2 - and the loop functions Mi, M 2 , M 3 and M 4 are given in Ref. [31]. The 
operators can be written as: 

u 1 — o 2 m s a L a s R r fJ.v , 

= ^rnJla»X P 4Gi ■ (52) 

One should also consider the C 7i5 7i9 contributions to 7^y, which are computed from 
Eq. (JH2|) as above described. 

The only dimension-four contribution relevant for this calculation is associated with 
the operator Qz, generated by the sdZ vertex which is mediated by chargino exchanges. 

G 6 

= ^— m 2 z tan 9 w Z ds s L ^^ d L + H.c. . (53) 
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The coefficient Z ds is given by 

Z ds = XtCo(x t ) + X t H (x qx ) , (54) 

where the first term refers to the SM Z° penguin contribution, evaluated in the 't Hooft- 
Feynman gauge, and the second one represents the SUSY contribution associated with 
At = (#23) l r(& '13) ] lr- Finally x t = mf/m^, x qx = m^/m x± and the loop functions are 
given by jSHj 



C n (x) = - 



x — 6 3x + 2 

7 H 9 In x 

x- 1 (x-1) 2 K J 



—x(x 3 — 6x 2 + 3x + 2 + 6x In x) 

Hn{x) = 



48(1 -x) 4 

(55) 

In USY scenarios, the SM prediction for e'/e, which receives contributions from the op- 
erators Oq and 8 , as well as from Z ds , is found to be 0(lCr 6 ). The supersymmetric 
contributions to e'/e could be dominated either by gluino or chargino exchanges. We 
decompose these contributions as follows: 



**{ 7 ) = Re ( 7 J + • < 56 > 

The first term is associated with gluino mediated diagrams and receives contributions from 
the "magnetic" operators in Eq. (|52|). The matrix element ((/T/t) 2 |O 7 |iC ) is proportional 
to the photon condensate, hence it must be very small. Accordingly, we shall neglect the 
C 7 contributions to e'/e. The matrix elements of the operator O g are given by 

((TTTT)JO n \K°) = 0, 



[fa) Q \O g \K°) = ^((n+n-)\O g \K°) 



3 1 11 m s p K _, 
2 167T 2 2 m s + rrid /J 



B q 1,2) , (57) 



where B^ 2 ^ ~ 1 and is the pion decay constant. Thus, the chromomagnetic contribu- 
tion to e'/e is 



]{vi - ] " 64^eWA WT^ d if ^ ml I,U 



(58) 



We now turn our attention to the chargino contributions. In addition to terms pro- 
portional to a single mass insertion (S^ll [3Tj . one has to take into account other con- 
tributions involving a double mass insertion, like those arising from the supersymmetric 
effective sdZ vertex. The term (e'/e)* can be separated into 



£ r\x f e'\* /e' N * 



Re ( — ) = Re ( - + Re I - . (59) 

£ J 1MIA V £ J 2MIA 
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As in Ref. j^I], the first term is given by 

(~Y = Im (j2 K *a2^1)LLK bl ) F £ /(x qx ) . (60) 

V E J IMIA V a,b J 

As before K denotes the CKM matrix, a, b are flavour indices and the function F e i(x qx ), 
with x qx = m| ± /m?, is given in Ref. |31j. This contribution is clearly dominated by 
(^12) ll- The second term on the r.h.s. of Eq. (|59|1 can be written as 



Re 



2MIA 



1.2 -J2, \rf\B { i /2) 



ImA t , (61) 



where the parameters i? s and are given in Ref. jSH]. 5f /2) is the non-perturbative 
parameter describing the hadronic matrix element and A t is given as 

A* = [(S^) LR (5^)U H (x qx ) , (62) 

where H has been defined in Eq. ()55|) . 

Now let us discuss in particular the role of each contribution within our model. As 
pointed out, the gluino contributions to e'/e can naturally saturate the experimental 
values if the trilinear couplings are non-universal, since the required value of the LR 
mass insertion is comparatively weaker, (lm(5f 2 )LR ~ 10~ 5 ). As mentioned above, it is 
possible to obtain such values and still saturate Sk, so the chromomagnetic operator O g 
can give a significant contribution to e'/e. Furthermore, and as aforementioned, in USY 
models the LL mass insertions are also enhanced, hence one might expect significant 
contributions from the chargino as well. 

In Fig. 0] we plot the gluino and chargino contributions to e' / e as function of the 
parameter 5 A that parametrises the non-universality of the trilinear couplings. 

5 A = (A 3i — Aii )2 i)/m = c — a (assuming a = b) . (63) 

The SUSY CP violating phases have been fixed as follows: ip a = (p c = and (pt, — 0.1. As 
in the previous figures, we have assumed tan/3 = 5 and m Q = mi/2 = 250 GeV. 

As can be seen from this figure, the chargino contributions with one mass insertion give 
the dominant contribution to e'/e. The chargino contributions with two mass insertions 
that arise from the SUSY effective sdZ vertex are also relevant and considerably larger 
than that of the gluino. 



4 Small CP phases versus large CP asymmetry aj/^Ks 

For a long time all experimental evidence of CP violation was confined to the kaon sector. 
However, recent results from the 5-factories have confirmed the existence of CP violation 
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5A 

Figure 4: Contributions to e'/e versus the non-universality parameter 8 A. The dotted 
line refers to gluino contributions, the dott-dashed and dashed lines to the chargino one 
(with one and two mass insertions, respectively). The solid line stands for the total 
supersymmetric contribution to e'/e. 



in the B sector as well. The fact that the asymmetry aj/^,K s is large may lead to the 
conclusion that CP cannot be viewed as an approximate symmetry. We will show next 
that this is not necessarily true. In fact, it might happen that models with small CP 
violating phases but with a considerably large flavour mixing can also account for the 
recent experimental results. As it was shown in previous sections, supersymmetric USY 
models are indeed an example of this class of models, where a large flavour mixing and 
small CP phases still succeed in saturating the measured values of the CPV observables. 
In this section we address some of the implications of supersymmetric USY models for 
the -B-system. 

Regarding the CP asymmetry of the Bd and Bd meson decay to the CP eigenstate 
J/ipKs, defined by 

aj/ ^ {t) = T(B d (t)^^K s ) + T(B d (t)^^K s ) = ^ Am ^ > ( 64 ) 

BaBar and Belle give the following results 

aj/^Ks = 0.741 ± 0.067 ± 0.034 (BaBar) , 

a J/i!Ks = 0.719 ± 0.074 ± 0.035 (Belle) , (65) 

which dominate the present world average, aj/^K s = 0.734 ± 0.054 [TU1 ITU I4H]. 
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In the case of the Standard Model, clj/^k s can be easily related to one of the inner 
angles of the unitarity triangle: 

V td V* 



V/iKs = sin 2/? ; /3 = arg --^ . (66) 



In the USY ansatz of Eq. (j!8j) . SM contributions to sin 2(3 are negligible, (9(10~ 3 ). In 
view of this, it is mandatory that new physics provides a large contribution to Bd — Bd- 
In the case of supersymmetric contributions to the AB = 2 transition, the off-diagonal 
element of the Bd mass matrix can be written as 

M 12 (B d ) = ^iy =2 l^) = M su [Bd) + M susy ^ (6?) 
2m Bd 

The supersymmetric effects are usually parametrised through 

rdC -Mf™( Bd y (68) 

where Am^ = 2|Mf 2 M (i? ( i)|r^. Thus, in the presence of SUSY contributions, the CP 
asymmetry parameter is now given by 

a jms = sin(2/3 + 26 d ) . (69) 

Hence, the measurement of aj/^Ks would not determine sin 2/3 but rather sin(2/3 + 26 d ). 

In many models, the SUSY contribution can be sufficient to account for the observed 
values of aj/^K s - However, this requires a large mixing in the LL and LR sectors, and 
this is indeed what we have in our USY example. As emphasized above, due to the large 
mixing in the USY textures, the non-universal A-terms have a non-negligible effect in 
the running of the squark masses and although we imposed universality of the scalar 
soft masses at the GUT scale, one finds considerable off diagonal contributions at the 
electroweak scale, particularly associated with mixtures of the third generation with the 
first two. Since in this class of models flavour mixing is not supressed, one finds that the 
(^13") ll and (5 13 u ) lr mass insertions (the most relevant in the aj/^Ks calculation) can be 
significantly large and saturate the experimental result as we show in detail below. In our 
numerical analysis, we took into account charged-Higgs, chargino and gluino contributions 
to Mf 2 USY (Srf), as given in Ref. [H]. The leading supersymmetric contribution to sin(2/3 + 
29 d ) is associated with the AB = 2 chargino mediated box diagrams. Instead of using 
the mass insertion approximation, as done in the previous sections, we choose to compute 
the full contibutions to M^ 2 (Bd)- The off-diagonal element of the Bd mass matrix is 
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accordingly given by: 



air 2 



16 3" 



2 

m . ± 

h,k=li,j=l Xj 



v a (r^ y„ x)Ll [rc r£ - v* 2 (r % ^ k), 6 



-tfcXj' X u h Xji X XiXj. 



} 



(70) 



where x a b = (m a /mb) 2 , /s d is the Bd meson decay constant, B Bi is the vacuum saturation 
parameter, and rj is the QCD correction factor. As before, V is one of the unitary 
matrices that diagonalise the chargino mass matrix and K is the CKM matrix. Yjj = 
diag(m u , m c , m t ) / (V2mw sin /?), H are the matrices relating the squark interaction and 
mass eigenstates, and the loop function G' is given in Ref. [33]. In Fig. El we present the 
correlation between the values of Ek and sin 29 d for tan (5 = 5 and mo = mi/2 = 250 GeV. 
The values of the parameters a, b, and c are randomly selected in the range [—3,3] and 
the phases are fixed as (p a = (p c = and <pi = 0.1. 
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Figure 5: Correlation between e K and sin 29 d . The crosses (+) correspond to values 
a = b= -3 and cG [-1.5,0.5]. 



From Fig. it is clear that one can have sin 29 d within the experimental range while 
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having a prediction for ek and e'/e compatible with the measured value. This result ap- 
pears as a characteristic of the kind of models under consideration. Despite the smallness 
of the phases introduced in these models via the Yukawa and trilinear couplings, they are 
sufficient to account for CP violation in the -ff-system, as well as the observed value of 
aj/^K s , which requires a phase of order one. 

It is worth mentioning that the non-universality of the A-terms is crucial for en- 
hancing the SUSY contributions. As emphasized above, for universal A-terms the mass 
insertions (S^ll and (S^ 3 )lr are, respectively, of order 10~ 5 and 10~ 6 , leading to a neg- 
ligible sin28d- When non-universal /1-terms are considered, we find that the chargino 
contribution displays a large phase, while the gluino and SM contributions are almost 
real and in general smaller than the chargino contribution. The Higgs and neutralino 
contributions are negligible in the cases under consideration. 

If we assume that the universality of the A-terms is broken by setting a = b 7^ c as in 
Eq. fl2Uj) . it is possible to find sizable values for the the SUSY contributions to aj/^K s - In 
Fig. El we considered particular values of a, b, c, namely a = b = — 3 and c G [—1.5, 0.5], 
and marked the associated contributions to sin 29 d and Ek with a cross (+). As seen from 
Fig. for this set of points, sin 29 d ranges from 0.99 to 0.3 while Ek takes values between 
(1 — 3) x 10~ 3 . Regarding values of the relevant mass insertions, the values of a, b, c 
corresponding to the crosses in Fig. we find that y] Im [(5£ L )f 3 ] |, a/| Im [(5^,) 13] | 
and a/| Im [(5£l)i3] Im [(5^)13] |, are of the same order of maginitude and range from 
0.08 to 0.2 as c varies from -1.5 to 0.5. The mass insertions Im [(52l)i3] I m [(^ll)^] |> 
VI Im [(5^)13] Im[(^ L ) 23 ]|, VI Im MlR I™ KKMl take values from 0.065 to 0.12 
in the same range. 

We notice that the values of these mass insertions are compatible with those given 
in Ref. j!2] to saturate the SUSY contributions to aj/^ Ks *. As a specific example, for 
a = b = —3 and c = —0.9 we find Ek = 2 x 10~ 3 , and sin 29 d = 0.7. The most relevant 
mass insertion for the gluino contribution is: (S 31 )ll = —0.0012 — O.O662, being the other 
contributions below 10~ 5 . The relevant mass insertions for the chargino contribution 
are (S^) LL = 0.0014 - 0.087i, (5'^) RL = 0.00051 - O.lli, {6% 2 ) LL = 0.011 + 0.047i and 
(5% 2 ) RL = -0.018 + 0.059i . 

Note that although we assume that A u = A d at the GUT scale, due to the effect of the 
large top Yukawa in the evolution to the electroweak scale, one obtains Im^^)^] « 0.1 
while Im[(53 1 ).Ri,] ~ 10~ 5 . Therefore, the chargino exchange gives the dominant contribu- 
tion to sin 29 d- 

To conclude this section, we address the AB = 1 process b — *> S7. As aforementioned, 

■"■The set of input parameters corresponding to Fig. [S] are (rriq) — 534 GeV, fi — 478 GeV, = 
210 GeV and the mass of the right-handed stop mj = 289 GeV. 
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USY models have large mixing, hence one might expect that the mass insertions {$23) ll 
and (S^lr receive large contributions, and this could lead to excessive values for the 
branching ratio (BR) of the b —>■ sj decay, above the experimental limit reported by the 
CLEO collaboration [33] : 

2.0 x 1(T 4 < BR (B — > X s7 ) < 4.5 x 10" 4 . (71) 

The supersymmetric contributions to the process are given by the one-loop magnetic 
dipole and chromomagnetic dipole penguin diagrams, which are mediated by charged 
Higgs boson, chargino, gluino and neutralino exchanges. 

It is known 03] that with universal A-terms, the gluino contribution is very small and 
the leading ones are those involving the charged Higgs and charginos. In this case, and 
for our choice of USY phases, we found that the total BR (B — > X s j) ~ 3.5 x 1CT 4 , 
which is essentially the SM central value. For the case of non-universal trilinear terms, we 
present just an illustrative example. Let us consider (for the structure given in Eq. PHjl ). 
a — — 1.7, b — c — —3, a point in parameter space associated with a nearly maximal value 
of sin(2/5 + 26 d ). In this case BR (B -> X s ^) ~ 2.8 x 10" 4 , still in agreeement with SM 
results, as discussed in Ref.jH]. 

5 Discussion and Conclusions 

In this work we have studied the implications of having universal strength of Yukawa 
couplings within the MSSM. This class of ansatze for the quark mass matrices is specially 
interesting since it provides a relationship between the observed pattern of quark masses 
and mixings and the possibility of having small CP violating phases. We have addressed 
the problem of CP violation in supersymmetric USY models, showing that the strength of 
CP violation due to the standard KM mechanism is too small to account for the observed 
CP violation, and essential contributions from supersymmetry are required. 

We have shown that the trilinear soft terms play a key role in embedding USY into 
SUSY. In fact, due to the large mixing and associated phases, the constraints from the 
EDMs on the SUSY parameter space are far more stringent than in the case of a stan- 
dard Yukawa parametrization. Although the effect of these phases is absent in the SM 
and in SUSY models with flavour conserving soft breaking terms, once we require a 
non-universality in the trilinear terms (which proves to be essential to saturate the CP 
observables) there are large contributions to the electric dipole moments of the neutron 
and mercury atom. We found that in order to satisfy the bound of the mercury EDM 
the A-terms should be matrix factorizable, and their phases constrained to be of order 
1(T 2 - KT 1 . 
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Within the region of the supersymmetric parameter space where one had compatibility 
with the EDM measurements, we have investigated the new contributions to both K and 
B system CP observables, and we have shown that in the present model it is possible to 
saturate the experimental values of Ek and e'/e. Gluino mediated boxes with LL mass 
insertions provide the leading contributions to £k, while e'/e is dominated by chargino 
loops, through LL flavour mixing. 

Finally, we considered the CP asymmetry of the B meson decay, a.j/^K s - It turns 
out that within this model the SM contribution to this asymmetry is negligible, while 
supersymmetry, namely through chargino exchange, provides the leading contributions, 
which are in agreement with the recent measurements at BaBar and Belle. 

In conclusion, we have presented an alternative scenario for CP violation, where the 
strength of CP violation originated from the SM is naturally small, due to the observed 
pattern of quark masses and mixing angles, which constrain all USY phases to be small. 
We have shown that in this framework new SUSY contributions are essential in order to 
generate the correct value of ex and e'/e, as well as the recently observed large value of 
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Appendix: USY and a local horizontal symmetry 

In this appendix we show how USY textures can be motivated by a horizontal symmetry. 
For instance, with the horizontal gauge group SU(3)h x U(1)h, the matter content of 
the MSSM can be assigned as {Q a ,L a } = (3, 1) and {u a , d a , e a } = (3, —1) , while the 
charges of the MSSM Higgs are given by {H ul Hc/\ = (1,0). The extra Higgs fields that 
may be used to break SU(3)h x U(1)h are <f) a = (3, 1) and <f) a = (3, —1), a = 1,2,3. We 
introduce an additional discrete symmetry, under which {H u ,Hd} — > {— H u , — Hd}, 
{<p a ,<f) a } — > {i<f) a ,i<f) a }, while all other matter superfields transform trivialy. This sym- 
metry prevents the presence of renormalizable terms involving quarks and leptons in the 
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superpotential. Therefore, the lowest dimensional SU(3)h x U{1)h invariant operators in 
the superpotential, which are responsible for generating the fermion masses, are given by 



W Yuk = KQ a u c b H u ^ + h d Q a d c b H d ^ + hiL a e c b H d ^-, (A2) 

where M is a scale much higher than the weak scale. The minimization of the scalar 
potential involving the fields (0 a , <p b ) can be derived in a similar way as in Refs. [2*T| I45j . 
Thus Y™ b = huj^r, where v a , v b are the vacuum expectation values of 0" = v a S ai and 
<f) b = v b 5bi (i = 1,2,3), respectively. In fact, for a particular choice of the soft terms 
associated to the fields (<j) a , (f) b ) one can assume that if v a and v b are given by v a = \v\e tlpa 
and v b = \v\e l,p b, then 

Yl = K^^\ with K = ^f. (A3) 

Similar expressions hold for Y d and Y l . The Yukawa couplings in Eq. (A3) clearly display 
the usual USY form, which was discussed in section |21 
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